Must Kohn-Sham Oscillator Strengths be Accurate at Threshold? 
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The exact ground-state Kohn-Sham(KS) potential for the helium atom is known from accurate 
wavefunction calculations of the ground-state density. The threshold for photoabsorption from this 
potential matches the physical system exactly. By carefully studying its absorption spectrum, we 
show the answer to the title question is no. To address this problem in detail, we generate a highly 
accurate simple fit of a two-electron spectrum near the threshold, and apply the method to both 
the experimental spectrum and that of the exact ground-state Kohn-Sham potential. 



I. INTRODUCTION 

Ground-state density- functional theory (DFT) [l|, [3, d, 
4] is enjoying more and more popularity for calculat- 
ing various atomic and molecular properties. The bal- 
ance between accuracy and calculation speed in DFT 
is achieved using an auxiliary Kohn-Sham(KS) system 
of non-interacting electrons. If the exact exchange- 
correlation energy were known as a functional of the den- 
sity, DFT would yield exact ground-state energies. 

In principle, all atomic and molecular properties are 
functionals of the ground-state density, including the 
properties of excited statesQ, but in practice only the 
ground-state energy functional has been usefully ap- 
proximated. The excited-state properties of the non- 
interacting Kohn-Sham reference system are often used 
to understand and even approximate those of the true 
interacting system, but in most cases this has no the- 
oretical justification. Thus the results of excited-state 
calculations with ground-state DFT must be carefully 
examined, since the KS orbitals and energies are (within 
ground-state DFT) artificial constructs designed only to 
reproduce the ground-state density. The more we under- 
stand about the differences between the KS system and 
the real system, the better we can determine whether an 
excited property of the KS system can be justified as an 
approximation to the real property. We study the exact- 
ness of the KS oscillator strength at the first ionization 
threshold in this paper. 

On the other hand, time-dependent density-functional 
theory(TDDFT) in principle gives several exact proper- 
ties of excited states 6] . Linear response TDDFT is a 
method that begins from ground-state DFT, and couples 
ground- state KS transitions to give the correct prop- 
erties of excited states. @, Q If we could use the exact 
time- dependent functional, the TDDFT method would 
exactly generate the properties of the real system from 
the results of the ground-state KS calculation of systems 
with non-interacting electrons. Thus our study of the 
exactness of the KS oscillator strength is converted to 
a question about the difference betweeen ground-state 
DFT and TDDFT. 

This may seem to be a simple problem, since ground- 



state DFT is not designed to give the correct oscilla- 
tor strength at the ionization threshold. The oscillator 
strengths can be extracted from strength of the poles 
of the linear response function, and the KS linear re- 
sponse function does not involve the Hartree-exchange- 
correlation(HXC) kernel(Refer to Eqn. {JTDJ and (JTTJ)), 
so there is no a priori reason to expect the exact KS 
system to give the correct oscillator strength at the ion- 
ization threshold. By 'exact KS' we mean the KS poten- 
tial as extracted from an extremely accurate ground-state 
density^, thereby avoiding the difficulty of distinguish- 
ing the effect of approximate ground-state XC function- 
als from that of KS-DFT itself. However, the ionization 
threshold of the exact KS system is equal to the ioniza- 
tion threshold of the real system, since Koopmans' theo- 
rem holds exactly for the exact KS system. [J] Thus this 
specific excited-state property is given exactly by the KS 
system, despite the lack of input from the Hartree and 
XC kernels. This is the only known direct link between 
real excited-state properties and their KS counterparts, 
and such links have proven invaluable in studying and 
understanding both ground-state and TDDFT jiol]. 

Given the usefulness of such links, and how both the 
strength and position of the the threshold occur at the 
same frequency, it is important to ask the title question, 
to see if some unknown exact condition might be lurking 
beneath the surface. To do this, we study one specific 
case. If, for He, we find definitively that the threshold 
oscillator strength is not given by the KS system, the 
answer is definitely no, and this cannot be true in general. 
If we did find it to match, the title question would remain 
open, and we would look for other cases and/or a proof 
of the equality. As we show below, the answer is indeed 
no. 

This has important consequences for the unknown ex- 
act XC kernel of TDDFT. To shift an ionization thresh- 
old, the kernel would need to be complex, with a branch 
cut at the position of the KS threshold. This is not the 
case for the first ionization threshold, but is for all higher 
ionizations. On the other hand, since we show that, typ- 
ically, the oscillator strength of the KS system is cor- 
rected by TDDFT, this means that, at the threshold, the 
Hartree- XC kernel must have some non-zero off-diagonal 



matrix elements. To understand this, showed that, 
in the absence of off-diagonal matrix elements, the KS 
oscillator strengths are unchanged by the action of the 
kernel. 



II. BACKGROUND 
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FIG. 1: KS and experimental single-electron oscillator 
strength of He near threshold. [U0 The ionization threshold 
is at 0.9036 Hartree. The bound-region spectrum is renormal- 
ized with the factor 1/n^ , where n / is the principal quantum 
number of the final state. 



In this paper we will show the answer to the title ques- 
tion is no. To prove that the KS oscillator strength 
doesn't have to be exact at the ionization threshold, we 
only need to provide a counter-example, so we study the 
helium oscillator strength spectrum in this paper. The 
helium atom is the simplest multi-electron system, and 
thus a theorists' favorite. Fig. [1] shows the photoab- 
sorption spectrum of Helium near the ionization thresh- 
old^. 6 eV), and the two curves are of the real helium 
and of the exact ground-state KS helium. Fig. [1] sug- 
gests the answer to the title question is no, but there 
could conceivably be near-degeneracies near the ioniza- 
tion threshold, and we wish to demonstrate that the os- 
cillator strength curve can be expected to be smooth near 
the ionization threshold explicitly. Hence we use a fit to 
explicitly show that the oscillator strength curves of the 
real helium and the KS helium are smooth across the ion- 
ization threshold, showing there are no near-degeneracies 
at the threshold. This allows the comparison of the value 
of the oscillator strengths at the ionization threshold of 
these two systems, showing that the difference between 
the spectra near the threshold of KS and real systems is 
inherent. Since our purpose is to understand the differ- 
ence between KS system and the real system, the fit is 
not done to the data points but to the general properties, 
such as the oscillator strength sum rules. We test our fit 
on the hydrogen oscillator strength spectrum, and then 
apply the fit to KS helium, real helium, and the result of 
using approximated TDDFT on the exact ground-state 
KS spectrum. 



In this section we provide a brief definition of notation 
and concepts used in this article. For any interacting 
electronic problem, the exact ground-state KS system is 
described by the KS equations 0: 



V 2 + i; s (r)Ui(r) = ei^(r), 



(1) 



where 



(2) 



Here and unless otherwise noted, we use atomic units 
so that energies are in Hartrees and distances in Bohr 
radii. The KS orbitals and eigenenergies are denoted 
as ipi(r) and and v s (r) is the KS potential, which 
can be separated into 3 pieces: external, Hartree, and 
exchange-correlation(XC) potential^, Q. The exact de- 
pendence of the XC contribution on the density is un- 
known and many approximation schemes are available, 
but in this article we use the exact value @. This is 
calculated by first obtaining the accurate density from 
a quantum Monte-Carlo calculation, then inserting the 
density into the KS equations and finding the potential 
that gives this density. Q 

The absorption spectrum in terms of photoabsorption 
cross section a is defined as below 1141: 



2tt 2 

Cr(w) = ^ fqt( W ~ U q) + O'contM, 



(3) 



where q denotes bound-to-bound transitions from state 
i to state f, f q is the oscillator strength of transition q, 
and <7 C ont which begins at uo = / is the spectrum of the 
continuum region. For bound-to-bound transitions, the 
oscillator strengths are defined as: 



(4) 



As defined in Eqn. ([3]), the spectrum comprises the 
discrete bound-to-bound transitions and the continuous 
bound-to-continuum transitions. We define a(oj) as the 
analytical continuation of cr con t(o;) for uj < I. This can be 
found easily by considering the bound oscillator strength 
as a continous function of uj q , yielding [3| 



a(co) = f{w)/ 



(5) 



In reverse, the usual oscillator strength for transition q 
Is — > np is given by 



dn) ~^ 



eh 



(6) 



where e is the energy of the np state. 
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The oscillator strengths are related to the dynamic po- 
larizability a(u>) by the following equations: 



a(u>) 



dr / d V zz'x(r, t';lj), 



= -^9f[a(o;)], 



(7) 



(8) 



where x( r , r '; is the linear response function of the real 
system, defined by the Fourier transform of the linear 
response function in time: 



X(r,t,r',t') 



5n[v ext ](rt) 



5v ext {r>t>) 



(9) 



«8xt [HO] 



The polarizability and the linear response function of the 
KS spectrum is defined similarly. The KS linear response 
function is related to the real(or 'exact TDDFT') linear 
response function by a Dyson- like equation [l5j]: 

X(r,r';cj) = x K s(r,r';w) + / d 3 ri J d 3 r 2 Xks (r, rj. ; w) 
x / Hx c(ri,r 2 ;w)x(r 2 ,r';w), (10) 

where Jhxc(i"j r'; uj) is the HXC kernel in frequency do- 
main, defined as the Fourier transform of the HXC kernel 
in time domain: 



/ H xc(r,r';t-t') = 



1 



5n(r',t') ' 



(11) 



The linear response function is also represented in 
Lehmann representation: [23j 



C(r,r';w) = lim ^ <^ 



ffoWsHM fl5Wffo(r') 



w — fi Q + %r\ uj + Q a + irj 



(12) 

where g a ir) = (^gs \n(r)\ and 17 Q = i£ Q — i? gs . 

Sum rules are moments of the oscillator-strength spec- 
trum, and they are related to various theoretical or ex- 
perimental physical properties of the ground-state atom. 
They are expressed with the following formula: [ltH [13] 



/■oc 

Sj = u){ cr s + / dui uPa(ui), 



(13) 



where s denotes the discreet Is — > np transitions and 
j is an integer. We only use — 2 < j < 2 in this article. 
These sum rules have simple relations to physical proper- 
ties, such as the ground-state density, polarizability, and 
kinetic energy, and thus they are easily calculated or de- 
termined from experiment. The specific relations we use 
are: 



S_ 2 = a(0), S_ 1 = -(|E i r j | 2 )o, S = N, 
Si = !(IE 3 Pjl 2 )o, S2 = Uzn(0), 



(14) 



where Z is the nuclear charge. Eqn. (fl3|) and Eqns. (|T4[) 
not only provide connections between the spectrum and 
several physical properties, but also imply that So, and 
S 2 are identical in the KS and the real spectrum, since 
the ground-state density in exact DFT is by definition 
equal to that of the real system. These equations also 
suggest the possibility of a fit which takes general physi- 
cal properties as input and is able to generate the entire 
spectrum for the H atom. 

As shown in Eqn. (|13p . sum rules of the bound region 
of the spectrum are calculated by a summation of the dis- 
crete Rydberg states, and there is no trivial formula for 
calculating the energies of these states in multi-electron 
atoms. In order to characterize these energies for the 
summation, we use quantum defect theory [14|. In quan- 
tum defect theory, the energy of the orbital with principle 
quantum number n in a multi-electron atom is expressed 
thus: 



EL 



1 



2(n - fi n ) 2 



(15) 



This expression is used in calculating the bound part of 
sum rules if the formula for /i is known. The quantum de- 
fect is a smooth function of energy, and can be very accu- 
rately approximated [18] by its Taylor expansion around 
fi = 0: 



,(p) 



(E)=J2^E\ E 



(16) 



i=0 



For Helium, this curve is essentially linear, so fj, ~ 
fj, + ViE, where /i = 0.0164 and m = 0.0289 for KS 
helium, and no = —0.0122 and /ii — —0.0227 for real he- 
lium. Inserting this expression into the E n formula and 
solving self-consistently yields highly accurate excitation 
energies [l9j 

To illustrate these features in an exactly soluable case, 
we use the hydrogen spectrum as an example in this ar- 
ticle. The exact form of the hydrogen oscillator strength 
is available using Eqn. (|4]):(n is the principal quantum 
number) 



0*ls — >kp 



256n ! 



(^)""/[3(n 2 -l) 4 ], (17) 



128 exp[/i(fc)] csch / [3(1 + k 2 



(18) 



where h(k) = { 7 r + 2tan- 1 [2k/ {k 2 - 1)] - 2ir6(k - 1)} /2, 
k = \f2e is the wavevector of the continuum wavefunc- 
tion, and 9 is the Heaviside step function. The 
wavefunctions of the bound states are energy-normalized 
by ri 3 / 2 , so the bound state wavefunctions and contin- 
uum wavefunctions agree with each other at ionization 
threshold. The hydrogen spectrum is shown in Fig. [5] 
We represent the bound transitions as by simple line 
segments whose height is d. 
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FIG. 2: Exact hydrogen spectrum. The ionization threshold is 
at 0.5 Hartree. Note that the bound-region spectrum has been 
renormalized so that it joins smoothly with the continuum- 
region spectrum. 



becomes exact when k — > oo. The transition dipole ma- 
trix element is evaluated at k — > oo limit. 



(ijjf \r cos#| ipi) 




h l—Z<p[{Q)\k~*/ 2 ,k 



The oscillator strength spectrum then decays as 



a(ui) 



2^2 



(21) 



(22) 



Eqn. (j2"2"|) implies the asymptotic decay of the oscil- 
lator strength only depend on the properties at the nu- 
cleus. For hydrogen and helium, Eqn. ([22]) is related to 
the electronic density by 



a(uj) 



8V2 



Z 2 n(0)w~ 



00. 



(23) 



III. HIGH-FREQUENCY LIMIT 

Fig. Q] suggests the KS oscillator strength and exact 
oscillator strength share the same asymptotic form. Real 
oscillator strength spectra of atoms decay as to' 
I20I |2l| . Here we derive the decay of the KS oscillator 
strength. 

The oscillator strength is related to the transition 
dipole matrix element |rcos#| by Eqn. (J4|. In 
the KS system, the matrix element is greatly simpli- 
fied, and can be written with one-electron KS orbitals as 
(ip{ |rcos0| ipi). For the absorption spectrum of the KS 
helium atom, the final orbital is a p orbital with wavevec- 
tor k, and the initial orbital is the Is orbital. In the 
high frequency limit of the absorption spectrum(w — > 00), 
k — ► 00 as well, and 4>)tp{r) is highly oscillatory, where (f> 
denotes radial wavefunctions. Then the matrix element is 
determined by the integrand near the nucleus. Thus the 
matrix element can be evaluated with the approximation 
of the initial KS orbital below: 



6i(r) = exp(— ar)[4>i(r) exp(cvr)] 

« exp(-ar) ( $(0)r + £[2a#(0) + #'(0)]r : 



(19) 

where <f>i is the spherical wavefunction of the initial KS 
orbital, and a is a positive real number characterizing 
the decay of the wavefunction. The cusp condition[22j 
holds in KS helium, so #'(0) = -2Z#(0), where Z = 2 
is the nucleus charge. Then is rewritten as 



4>i(r) « exp(-ar){r + (a - Z)r 2 }<j)[{0). (20) 

In k — > 00 limit, only the —2/r Coulomb well in the KS 
potential is important to (/>k p . Then </>k p is approximated 
with hydrogenic wavefunctions, and the approximation 



For hydrogen, the coefficient of the u> 7 / 2 term is 
8%/2/3tt. Eqn. ^ and ^ give the correct result. With 
these equations, the asymptotic behavior of the oscillator 
strength spectrum is determined. The discussion of the 
high-frequency part of the KS helium oscillator strength 
spectrum can be extended to other KS atoms easily, as 
the KS system is an one-electron picture. Following sim- 
ilar procedure as described here, it can be easily verified 
that the high-frequency part of the KS oscillator strength 
of other KS atoms can also be expressed in terms of the 
density at the nucleus. As only s orbitals has non-zero 
contribution to the density at the nucleus, one would ex- 
pect that Eqn. (|2"3"|) also holds for other atoms, using the 
corresponding Z and n(0). 
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FIG. 3: KS and experimental single-electron oscillator 
strength and fit curve of He near threshold. [l2l [l3l| The ioniza- 
tion threshold is at 0.9036 Hartree. The curves are converted 
from g(x) fit(Eqn. I|27[l). The upper curve represents the 
exact KS Helium oscillator strength data and fit curve, and 
the lower curve represents the experimental Helium oscillator 
strength data and fit curve. 

The half-power decay of Eq. (23) differs noticeably 
from the decay discussed by van Leeuwen, but here we are 



4 



considering the imaginary part of the response function, 
whereas there it is the real part of the response function. 
We are currently investigating relationship between the 
two in the general case. 



IV. FITTING OF 1- AND 2-ELECTRONS 
SPECTRA USING SUM RULES 

We fit the oscillator strength spectra to answer the 
title question. Since we want to study the near-threshold 
behavior of the oscillator strength spectrum, the position 
of the ionization threshold is treated explicitly in our fit. 
We define x and g(x) as 



/(*) 



(24) 



The fit has to satisfy a few criteria to generate the cor- 
rect shape for the oscillator strength spectrum. The fit 
is employed to study the exactness of the KS oscillator 
strength at the ionization threshold, where the fit needs 
to have the correct series expansion. We take the expan- 
sion of the hydrogen oscillator strength at the ionization 
threshold: 



a(u^I) = c + d(cj - I) + c 2 ((J - I) 2 + 



(25) 



We assume the fit formula has the same expansion near 
the ionization threshold. This assumption is justified by 
the following consideration. Near the ionization thresh- 
old, the oscillator strength spectrum of real helium is 
determined by the Rydberg states, which resembles the 
hydrogenic states. The KS helium is a system with non- 
interacting electrons, so the oscillator strength spectrum 
resembles that of one-electron systems. 

We use the fit to show that the oscillator strength spec- 
trum around the ionization threshold is smooth, that no 
near- degeneracies exist around the ionization threshold. 
Thus the fit also need to accurately generate the entire 
oscillator strength spectrum, including both discrete and 
continuum regions, so that the conclusions from the fit 
are convincing. To generate the correct continuum spec- 
trum, the fit need to have the correct series expansion 
when io — ► oo. As in Sect. IIIIl the asymptotic series 
expansion of helium has the same form as hydrogen: 



ct(w 



d\ui 



i 2 ui 



d 3 oj~ 



(26) 



The shape of the g(x) function is shown in Fig. [5] Our 
g-fit formula is 



g(x) = a + b[l — exp(— cx)] + dyje + x, 



(27) 



where a, b, c, d, e are fit parameters. Note that aside 
from giving the correct series expansion at the ioniza- 
tion threshold and asymptotically, the form does not have 
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FIG. 4: KS and experimental single-electron oscillator 
strength and fit curve of He.[l^. [13 ] This figure show the 
overall shape of the oscillator strength curves. The solid dots 
and curve represents the exact KS Helium oscillator strength 
data and fit curve, and the cross dots and dashed curve rep- 
resents the experimental Helium oscillator strength data and 
fit curve. 
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FIG. 5: g(x) of Hydrogen(Eqn. ([Ml)). The ionization thresh- 
old is at x = 0. 



other explicit physical motivation. It is solely designed 
to recover the shape of the oscillator strength curves. We 
determine the parameters by the process below. 

The important points are fixed on the fit. Since we use 
the fit to study the oscillator strength around the ioniza- 
tion threshold, we fix the value and the first derivative of 
the oscillator strength at the ionization threshold. The 
asymptotic coefficient in Sect. IIIIl is not fixed, but used 
as the initial point of search. The remaining three param- 
eters are determined by applying oscillator strength sum 
rules(Eqn. (fT5)) ) to the fit curve. We evaluate the sum 
rule of a fit curve by adding the contributions from the 
discrete transitions and that from the continuum. For 
the discrete region, we calculate the frequency of a tran- 
sition with quantum defect theory(refer to Eqn. (fl5l) 
and (fl6)) h The oscillator strength of the transition is 
then evaluated with the fit formula(with a certain ini- 
tial choice of parameters). We add the contribution of 
different discrete transitions up to n — 1000. For the 
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continuum region, we carry out a numerical integration 
over the entire continuum. 

The exact values of the oscillator strength sum rules 
are available for both KS and real helium, since these 
sums are related to various physical properties (refer to 
Eqn. I|14p). To fit the oscillator strength spectrum, we 
choose an initial set of the parameters. Only three fit pa- 
rameters are independent, so we choose three sum rules 
to fit. We minimize the difference between the sums eval- 
uated on the fit curves and the exact sums obtained from 
physical properties by varying the three parameters nu- 
merically. The search ends when the accuracy of the fit- 
ted sums reach a predetermined goal. In our application, 
the difference between the sums of the fit and the exact 
sums is smaller than 10~ 8 . The accuracy of the fit is also 
checked by evaluating the unused sum rules (Table. HJ). 

The fit can use two to four sum rules depending on 
how many points are fixed in the beginning. Applying 
more sum rules increases the overall accuracy of the fit, 
but the process of numerically fitting sum rules becomes 
more difficult. All results in this paper are obtained with 
three sum rules. With Eqn. (j2"Tj) . the sum rules of the fit 
curve can be written out in terms of the parameters. 



i 

0.8 
0.6 
0.4 
0.2 






ALDAfit 

ALDA data 
Experiment data • 
1st ionization threshold 



20 



30 



40 



50 



co(eV) 



FIG. 6: Exact /ALDA oscillator strength and fit curve of He. 
These oscillator strength data are obtained from an ALDA 
calculation with exact KS ground state. We use a box code[l9l] 
to calculate these data. There is a kink in our data near the 
ionization threshold, because the continuum near the ioniza- 
tion threshold mixes with higher Rydberg states, which are 
not well-described by the box code. 



TABLE I: Sum rules from ct(oj) and g[x) fit 



+I j - 3 [a + b + b{j-3) exp (2c/) Kt-j (2c/)]/(3 - j)} 



s ? s = 3 E £ 3 { a + b t 1 - ex p W\ + d V^W} ■ 



(28) 







S-2 




Si 


s 2 


H 
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He KS a 


g-fit 0.7563 
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g-fit 
exact 
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g-fit 
exact 


0.6912 c 
0.6912[26] 
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0.7957 


2.0414 
1.9114(27] 
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where j3 = (-^ - Hi/n 2 + n) , 7 = (/ - /3 2 /2)~ 
Ho and are the parameters in the quantum defect fo- 
rumla(Eqn. (I16p ). B is the incomplete beta function, and 
E is the exponential integral function [23], 

With Eqn. (|2T)l. we obtain the oscillator strength 
curves of KS Helium and real Helium. We also apply 
our method to the ALDA Helium(with exact KS ground 
state) as the first step of studying the threshold behav- 
ior in TDDFT (Fig. [5]). The comparison of results and 
figures of oscillator strength curves are shown in Sect. IVl 
and in Fig. [3] 

Note that the fit is not designed to be used as an inter- 
pretation tool, but to recover the shape of the oscillator 
strength spectrum. Thus comparing the fit parameters 
of different curves(exact KS, ALDA, and experimental) 
is largely meaningless as there is no visible trend. An ex- 
ception is the fit parameter d, which describes the shape 
of the asymptotic part of the oscillator strength curve, as 
it is related to the coefficient of the leading term(w~ 7 / 2 ) 
of the asymptotic expansion of the oscillator strength. 
The fit parameters of related systems are provided in the 
supplementary material. pHj 



a All the sums are converted to corresponding single- 
electron sums 

'The expected value of experimental data are listed in [171 
c This sum rule is a constraint. So is always a 
constraint. (So = 1 for all systems after converted to single- 
electron model) 

d The expected value of S— 1 and S2 are calculated from the 
exact helium density [j| 



V. RESULTS 

The g-fit curves of the KS helium and the real helium 
are shown in Fig. [3] and [5] The fit is very accurate 
in the entire range of u>. The accuracy is also checked 
with the unused sum rules, listed in Table. U The re- 
sults of the hydrogen atom are listed as a reference, and 
it shows that the inherent error of the method is small. 
With these curves, we explicitly show that the oscillator 
strength spectrum is a smooth curve around the ioniza- 
tion threshold, and thus the oscillator strength of the 
exact KS helium is not that of the real helium. The au- 
toionizing resonances in real helium are not included in 
our fit. but the fit is still accurate even near the reso- 



G 



nances (Fig. [7]). The errors in the sum rules are small, 
so the fit curve can be used as a background for studying 
these autoionizing resonances, and the pure resonance 
peaks can be obtained by subtracting the fit curve from 
the experimental spectrum. 

One reason for the good performance near the res- 
onances is that the autoionization resonances occur at 
relatively high frequencies, so their contributions to the 
smaller sum rules are neglectable. The other reason is the 
shape of the autoionization resonances in He is asym- 
metric, which have both a dip and a peak in the reso- 
nance region. [28, 29] The contribution of these two parts 
to the sum rules cancels, so the values of the sum rules 
are not influenced by the autoionization resonances too 
much(even for S2), and thus the fit accurately generates 
the oscillator strength curves for He. 
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FIG. 7: gf-fit of the experimental helium oscillator strengths 
near the autoionizing resonances. 



VI. CONCLUSIONS 

If the answer to the title question had been yes, then it 
would yield a strong exact condition on the XC kernel in 



TDDFT, which many approximations would fail. Thus 
we studied the title problem. We have shown that Kohn- 
Sham oscillator strength of He is not exact at the ioniza- 
tion threshold (even though the position of the threshold 
is exact), and so the answer to the title question is no. 
This implies that the Hartree-XC kernel in TDDFT has 
non-zero off-diagonal matrix elements at the threshold, 
and simple approximations such as the single-pole ap- 
proximation are insufficient in this region. 

We also developed a numerical fit to generate the spec- 
trum near the ionization threshold from a few physical 
conditions such as sum rules. The fit is accurate for 
all frequencies due to the smoothness of the oscillator 
strength near the threshold, but also works well for the 
spectrum far from the ionization threshold due to the 
correct asymptotic behavior. The fit is not physically 
motivated, but is a simple accurate representation of the 
curves. 

These results are not general since we only studied 
atoms with one or two electrons, and multi-electron res- 
onances are ignored as in Fig. [3j However, obvious 
generalizations can be performed for atoms with more 
electrons since we only use the general properties (the 
asymptotic behavior, the value and first derivative of the 
spectrum at ionization threshold, and sum rules) in our 
method. Thus multi-electron resonances can be dealt 
with by subtracting their contribution from sum rules, 
and thus our method can be extended to other atoms by 
following the methods of Sect. IIIII and [TVl 
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